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,.0 ■ Abstract 

We study thermodynamics of an ideal gas in Doubly Special Relativity. New type of special functions (which we 
call Incomplete Modified Bessel functions) emerge. We obtain a series solution for the partition function and derive 
thermodynamic quantities. We observe that DSR thermodynamics is non-perturbative in the SR and massless 
limits. A stiffer equation of state is found. 
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O" 1 Introduction 

Vh ■ 

£J) Attempts to combine gravity with quantum mechanics in search of the theory for quantum gravity always seem to 

give rise to the Planck length (jp = \J~^jr \ that provides the scale at which the quantum effects of gravity will show 

\ up [1-11]. The existence of such a length scale is in conflict with the equivalence principle because observers in 
different inertial frames will not agree on lp due to the Lorentz-Fitzgerald contraction. It has been shown that it is 
possible to still have equivalence principle by deforming Special Relativity (SR). These classes of theories fall under the 
qq ■ name Doubly Special Relativity (DSR) [12-15]. In DSR, apart from the constancy of speed-of-light scale, the Planck 
\ length lp or equivalently Planck energy k is also constant under coordinate transformation from one inertial frame 
qq ■ to another. This leads to modification in the dispersion relation. Consequences of the modified dispersion relations 
f^*) | on the thermodynamics are being studied extensively to infer the effect of Planck scale physics [16-22]. The effect of 
~l ■ modified dispersion relations in loop-quantum-gravity on black hole thermodynamics was studied in [16]. The same 
as a Lorentz violating phenomena on the thermodynamics of macroscopic systems (like white dwarfs) [17-19] and as 
J> ■ a noncommutative phenomena on cosmology and astrophysical systems [21,22] have also been studied. Moreover, 
photon gas thermodynamics in the context of modified dispersion relations [20] and DSR [23] are being investigated. 
In [23] the effect comes solely because of the presence of a maximum energy scale as the photon dispersion relation 
C$ \ remains unmodified. 

The present paper aims to study the thermodynamics of an ideal gas consisting of massive particles in DSR scenario. 
Both the modification in the dispersion relation of the constituent particles and the presence of a maximum energy 
scale are expected to contribute to new effects. DSR transformations can be of several type. In this paper, we follow 
the formulation of [24, 25] where the modified dispersion relation becomes 



- r— I 

X 
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e 2 - p A =m A [I - -) . (1.1) 

As < (1 - |) < 1, the energy of a particle with a given momentum decreases in DSR. This has consequence on the 
thermodynamics as we will see in § 3. The parameter m can be called "invariant mass" as it remains invariant under 
a DSR transformation. Note that in contrary to the SR case, m is no more the rest mass energy of the particle. To 
get the rest mass energy mo, we put p = in (1.1). We get two expressions for mo, namely 

±m 
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The two solutions are connected by the redefinition of the parameter m — > — m. Henceforth, without any loss of 
generality we use 

The physical world is characterized by E < k [25]. In this sub-Planck regime (E p= q = too < k), the positivity of rest 
mass (mo > 0) restricts the range of the invariant mass to < m < oo. Thus, in (1.1), we have < p, E < n and 
< to < oo. 

We study the thermodynamics of an ideal gas in DSR setup. We obtain a series solution for the partition function 
and compute the various thermodynamic quantities. We show that our results go to the standard results in the SR 
limit (k — > oo) [26] as well as in the massless DSR limit [23]. 

2 The Partition Function 

We consider a gaseous system of non-interacting particles obeying Maxwell-Boltzmann statistics whose macrostate is 
denoted by (N, V, T) where N is the number of particles in the system confined in volume V at a temperature T. In 
the canonical ensemble the thermodynamics of this system is derived from its partition function [27] 

Z N (V,T) = Y,eM-PE], (2-1) 

E 

where /3 = and ^ denotes sum over all the energy eigenvalues of the system. The total energy E of the system 

E 

can be written in terms of single particle energy e 

E = J2 n ee, (2.2) 

e 

where n e is the number of particles in the single-particle energy state e and satisfy the following condition 

5> e - N. (2.3) 

We can rewrite as 

Z N (V, T)=J2 SRI exp[-/3 J2 (2-4) 

{nj e 

where g{n e } is the statistical weight factor appropriate to the distribution set {n E }. The summation goes over all 
distribution sets that conform to the above restrictive condition (2.3). For Maxwell-Boltzmann statistics, it can be 
shown 

Z N (V,T) = ± i [Z 1 (V,T)] N , (2.5) 
where Z\ is the single particle partition function given by 

Z 1 (V,T) = J2^Pi-M- (2-6) 

e 

While for ordinary spacetime, it is easy to show that in the large volume limit one can replace the sum by an integral 

E -* p / d ^ ( 2 - 7 ) 

for more exotic spacetimes the measure of integration is expected to get modified 

d 3 p -+ f(p)d 3 p. (2.8) 
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Hence putting together (2.6), (2.7), (2.8) and taking h = ks = 1 we get 

V 



Zi (V,T) 



(27T)' 



■ r d 3 P f(p) 

J P =o 



exp[~P (s - mo)]. 



(2.9) 



Note that in accordance with standard practice, we have subtracted the rest mass Too from the relativistic energy e 
of the particle. Although there have been few attempts [19,28,29], the form of f(p) is far from settled. Assuming 
isotropy of spacetimc we may take f(p) = f(p). For a possible deformation of the integration measure, f(j>) should 
be expandable in Taylor series in ^ 



with ao = 1 since in the limit k — > oo, f(p) — > 1. Hence Zi(V,T) becomes 

T /" Aft OO 

(2tt) J p= o n=0 v « 7 

= --I 

where Z\ (V, T) is the single particle partition function with unmodified measure 

V 



Z°(V,T) 



(2n)- 



/ d 3 pexp[- 
J P =o 



P(e- too)]. 



(2.10) 

(2.11) 
(2.12) 

(2.13) 



The derivation of (2.12) from (2.11) involves two steps. Firstly, the interchange of the summation and integration 
which is allowed if (see theorem 1.38 of [31]) 

£^i> E "" pM(£ -" , ° )l - £^(--|)" z .°< i ' r ><°°- < 2 -»> 

Secondly, writing / d 3 p s n exp[— (3 (e — mo)] as (mo — f p — ex P[ — P ( £ ~ too)] since the former can be 

written as 4*7r f* =Q dp p 2 e n exp[— p (e — m )] whose integrand p 2 e n exp[— p (e — mo)] remains to be continuous and 
bounded for p e [0, k],P e [0, oo] (see § 5.12 of [32]). Hence our problem has boiled down to solving the integral in 
(2.13) where e and p are related by the modified DSR dispersion relation given in (1.1). The solution of Z° in the 
massless case obtained in [23] is 



7° - 



2V 



(2tt) 2 /3 3 



(2-e' fSK {p 2 n 2 + 2pn + 2)) 



(2.15) 



The term with e"^ K makes Z\ ml non-analytic at - = 0. We anticipate that even when too 7^ 0, Z\ continues to 
be non-analytic at - = and hence does not admit a straightforward Taylor series expansion in -. Thus in order 
to find the leading order deviation of DSR thermodynamics from the SR case, one would require a non-trivial series 
expansion. 

2.1 Solution of Z\ 

Changing the variable from p to e in (2.13) we get 



z°Av,P) 



2V 
(2^7 



exp[/3m ] 



f e+-(l--) e 2 -m 2 (l~ £ - 



1 1/2 



exp[— ps]de. 



(2.16) 



We now consider three different regions of values of to: 
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2.1.1 Case I: m = k 

In this case the partition function reduces to 

Z°{V,T) 21 



2 k 3 / 2 exp 



(27T) 

2\/2F 



/3k 



3/2 



A/2 
3 /3/t 

2'~2~ 



de(2e- k) 1/2 exp [-/3e] 



where j{a,x) is the Incomplete Gamma Function (see (6.5.2) of [30]). The factor s + ^- (l — ^ ) 
this simplifies the integral in (2.16) yielding a simple analytical form for Z\(y, T). 



(2.17) 

reduces to n and 



2.1.2 Case II: k < m < oo 

v 2 



We put t 



(t) - X J - f in ( 2J6 ) to § et 
2Vm 3 



z\ = - 



(2nf [(f) 2 1 
2Vm 3 

(2tt) 2 



3/2 



3/2 



exp 



exp 



(3m - 
f3m Q - 



i'm 2 



/ — K/m 
dt t(l -i 2 ) 1/2 cxp[-/3'mt] 

V m/J 



(2.18) 



v - 



I*(x,y)= I dt t(l-t 2 ) 1/2 exp[-a;t]. 
J-y 

We define Incomplete Modified Bessel function I v (z,y) of order v 

I v (z,y) = / —jr (|)" [ (l-t 2 y-^cxp[~zt]dt [Rev>0,\argz\<^], 
such that I*(x,y) = -^[f 7i(a;, j/)]. In particular for y = 1, using (3.387 (1)) of [33] and (9.6.26) of [30] we get 



(2.19) 



(2.20) 



I* (3,1) 



■nh{x) 



(2.21) 



where -^(a:) is the 2nd order Modified Bessel function. In the limit m — > 00 (m — > k) one gets from (1.1) 1 



e 2 -- 2 



£- = (l--) Vpe [0,/s]. 



Thus the total energy E of the system becomes E = Nk and the thermodynamics simplifies. Entropy can be computed 
by counting the total number of microstates Qn available to the system 



n N 1 

N M Nl 



K d*p 
p=o ^ 3 



1 / 2Vn J 



N 



NI\3(2tt) 2 J ' 

where fli is the total number of microstates available for a single particle. Thus the entropy S of the system is 



(2.22) 



5 = In 



1 / 2Vk 3 



N' 



Nl \3(2ttY 



(2.23) 



1 We would like to thank Diptiman Sen for pointing out this interesting case. 
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The first law of thermodynamics in this case becomes 

dE = -PdV + ndN. 



(2.24) 



Note that the usual term TdS has been dropped as from (2.23) it is evident that S is a function of N and V alone. 
The pressure of the system is zero asP=- w|/v = ^ wm l° the chemical potential is fi = §w\ v = K - 



Equation (2.13) can now be easily integrated to give 



2V k 3 
(2^fT : 



which gives the limiting behaviour of I* (/3'm, ^) using (2.19) and (9.6.7) of [30] 

I* 



13m K \ m ^oo 1 



2.1.3 Case III: < m < k 



We put t = j- 



1 - (f ) J + T in ( 2 - 16 ) to § ct 
2Vm 3 



2Vm 3 



1 - (2i 



3/2 



3/2 



exp 



cxp 



/3m 



/3m 



B"m 2 



/K/m 
dt t(t 2 -l) 1 ' 2 exp[-p"mt] 

AT* (/3"m, 1) - K* ((3"m,- 



where B" = — and 
l-(f-) 

/>OC 

K*(x,y) = dt t{t 2 - l) 1 / 2 exp [-art]. 
As in Case II, we define Incomplete Modified Bessel function K v (z,y) of order f 

K v (z,y)= r( J + 1) (-) y (t 2 -l)'^^exphzt]dt [ife i/ > --, |ar ff z| < -], 
such that K*(x,y) = — ]• In particular for y = 1, using (9.6.23) and (9.6.26) of [30] we get 

*•(*,!) = ^>, 

where if 2 (a;) is the 2nd order Modified Bessel function. 

We shall now obtain the series solution of K*(x 7 y). We rewrite (2.28) as 



K* (x 



,y)= r 

•>y 



1/2 



(2.25) 



(2.26) 



(2.27) 



(2.28) 



(2.29) 



(2.30) 



Inside the integral t > y and for y > 1 (which is a valid assumption for the case of our interest) the factor (l 
can be expanded in series of -k to get 



K* (x 



,y) = / 

Jy 



dm 



with 



d Mt = t 2 e~ xt dt 



(2.31) 

1 \V2 

(2.32) 
(2.33) 
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and 



where 



fr(t) — fir ' 

(Q-|)(l-i)...(r-l-l) 
r! 



(2r-2)! 
'22»-ir!(r-l)!' 



(2.34) 
(2.35) 



°° /-OC 

Now the integral and the summation in (2.32) can be interchanged if / dyU t |/ r (t)| is finite (see theorem 1.38 

of [31]). Now as t r is —we for all r > 1 we have \f r {t)\ = —fr(t)- This allows us to interchange the summation and the 
integral if the final series is converging. So we get 



K* (x, y) = M - -Mi + ^ i r M r , 



r=2 



if the above is a converging series (see Appendix A for convergence of K* (x,y)). Here 

/■OO 

M r = / dt t^-^e"*' 
■'1/ 



(2.36) 



(2.37) 



for r = 2,3,.... 

M and Mi can be easily calculated to be 



M = CXp( 3 Xy) ((xyf + 2xy + 2), 



Mi = 



exp (-xy) 



Now, changing the variable to t' = xt in M r for r > 2 we get 



M r = a; 2r - 3 / dt'- 



/\2r-2 ' 



(*') 



Taking e * as first function, if we do the integration by parts again and again we finally get 



(2.38) 
(2.39) 

(2.40) 



„2r-3 



2r-3 



M r . 



(2r - 3)! 



(-*) 



fc-i 



(2r-3)(2r-4)...(2r-2-fe) Vz// 



2r-2-fe 



(2.41) 



Here E\{x) is the Exponential Integral (see (5.1.1) of [30]). A similar attempt to obtain the series solution of I*(x, y) 
fails. 

Although we obtain the solutions of Z\ in three different regions of values of to, Z\ can be shown to be smooth in 
to (see Appendix B) and hence we do not expect any phase transition like thermodynamic discontinuity as we vary 
to. We use the continuity of Z\ to obtain the limiting behaviour of I*(f3'm, — ) and K*(f3'm, — ) as m — > k. From 
(2.17), (2.18), (2.21) and (9.7.1) of [30] we obtain the leading order behaviour of I* (/?'ra, ^) asTn -> k+ to be 



(2.42) 



with (5 = f -1. For to — ^ k~ , using (2.17), (2.27), (2.30) and (9.7.2) of [30] the leading order behaviour of K* (j3"m, ^) 
turns out to be 



K*(p"m,-) -A 

TO (/3k) j /' 



-e 3 / 2 e- 



I ■) , _ / 3 fin 



(2.43) 



where e = 1 — — . 
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2.2 Leading order deviations 

Having obtained the series solution of Z\ in Case III, we shall now obtain the leading order corrections from the 
massless and the SR cases. 

2.2.1 Leading order deviation from the massless case 

Thermodynamics of a photon gas in DSR with dispersion relation (1.1) and unmodified measure has been worked out 
in [23] . Here we calculate the deviation of single particle partition function from that of a photon gas. On expanding 
Z\ in rj = ^ with mo — > (assuming k to be finite) and using (9.6.10) and (9.6.11) of [30], we get 

Z l = Z lml + Z lmlcom (2.44) 

where Z® ml is the single particle partition function of photon gas in DSR scenario with unmodified measure [23] and 

= ^^(2-e^»(/^ + 2/fo + 2)), 



2V (Bk) 4 

Zlmlcorr = ( ^ )2/?3 V Hi) (v* + Otf)) + (f} K Z° lml ) V + Otf). (2.45) 



Note that the correction due to mass of the constituent particle is non-perturbative in nature as the first term in 
z imicorr which contains ln(rj) is the non-analytic piece and does not allow a trivial Taylor series expansion at rj = 0. 
We can rewrite (2.12) as 

Zl = Zi m i + Z lmlcorr, (2.46) 



where 



n=0 



Z lmlcorr = ZUorr + ^r^l^^-Vn.^^). (2.47) 

n— 1 \ y y / 

The above leading order behaviours have been plotted in Fig 1. For our choice of parameters they match with the 
numerical plots up to ^ — 0.012. 

2.2.2 Leading order deviation from the SR case 

On expanding Z\ in r\ = ^ with n — > oo (assuming m to be finite) , we get 

Z l = Z 1SR + Z \SRcom (2.48) 

where Z\ SR is the single particle partition function in SR and Z\ SRcorr is O (rj). 



_ 2Vm% 0mn K 2 {0m o ) 
2Vml el 3 " 10 e - /3m °/ T ' 



(2.49) 



Note that the DSR correction is non-perturbative in nature as the first term in Z® SRcorr which contains e /3m °/ 1 ' is 
the non-analytic piece and docs not allow a Taylor series expansion at rj = 0. This is a novel feature in DSR as we 
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Figure 1: The single particle partition function is plotted vs mo (left) and - (right). The plots for the unmodified 
measure are in solid lines while those for modified measure are in dashed lines. The plots in massless (left) and SR 
(right) cases are shown in red. The leading order behaviours are plotted in blue. The numerical plots for (2.27) and 
(2.12) are shown in green (solid and dashed respectively) for comparison. Different values for the parameters in natural 
units are as follows: V = 10 35 ,7V = 10 25 ,T = 0.01, a = l,a t = 0.2, a 2 = a 3 = ... = and k = 1 (left), m = 0.01 
(right). 



know that SR thermodynamics is perturbative in the non-relativistic limit: 

7° 4nV ( fc g TTO 0) 3/2 -l/2„«r^ („\ 

Z 1SR = — u e K 2 {u) 



/ 2nmok B T\ 3/2 ( 15 ,„ / 7 l^ 2 



-*"( 1 + s + ((i)j)- (2 ' 50) 

where u = and Zimr is the single particle partition function in the non-relativistic case. We can rewrite (2.12) 

Z\ = Zisr + ZiSRcorr, (2-51) 



k B T 

as 



where 

Zisr = Z^ SR , 

7 _ 7 , a l m yQ «1 dZ ° SR 

^lSRcorr — ^lSRcorr ^ ~ — ^lSR ~ go — • \4.Q£) 

The above leading order behaviours have been plotted in Fig 1. For our choice of parameters they match with the 
numerical plots up to ^ ~ 0.08. Having obtained the leading order correction to Z\ due to DSR, we shall now 
compute its effect on the various thermodynamic quantities. 
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3 Thermodynamic quantities 



The free energy F, pressure P, entropy S, internal energy U, internal energy density p and heat capacity Cy are 
defined as 

'Z£ 
N 



F 
P 



-\\a{Z N {V,p,m )) = ~N{ln 



fdF 
[dV 



1 



(3.1) 
(3.2) 



N,T 







s = 


-(-) 




\dTJv,N 


u = 


F + TS 




U 


p = 


V 


C v = 




(-) 



(3.3) 
(3.4) 
(3.5) 

(3.6) 

The above quantities for Case I can be found by using (2.17) and they have been plotted in Fig 2. Now we shall obtain 
the leading order thermodynamics for Case III. If we denote the free energy, pressure, entropy, internal energy, internal 
energy density and heat capacity obtained in the SR or massless cases by F , Po, So, Uq, po and Cyo respectively, and 
write Z\ = Z w + Z\ corr , where [Z w , Z\ corr ) = (Zisr, ZisRcorr) or {Z\ m i, Z lrn i corr ), we have 

2^ 

F - — - ~~ + O 



F 
S 



N 

F - —In [ 1 



J lcorr 



Zw 



N Zi cori 
J Zio 



'Z u 



Zio 



(3.7) 



N 



PN dZ 1 



Z l dp 
NZi corr PN ( dZi corr 
Z\o 



U 



C v = 



P = 




(3.8) 
(3.9) 



n d 2 z lc 



no 



d 2 p 



(3.10) 
(3.11) 



where n — ^ is the number density. The correction to F depends on the ratio of Z\ corr and Z w and is independent of 
the volume V of the system. Note that this is true to all orders. Hence, the pressure P of the system which is defined 
as P ~ — (ffOjvT S e * s no correc ti° n: 



P = Po- 

The equation of state in SR is (see (8.128) and (8.134) of [26]) 

Psr 



P S R 



which gives the following DSR equation of state 



P 



ZlSRcorr p _|_ n S^lSHco 



ZlSR 



O 



Z\SR,corr \ 
. Z 1SR ) 



(3.12) 
(3.13) 

(3.14) 



9 



1.4X10 - 



1.2 x 10 26 - 




Figure 2: The plots for the unmodified measure are in solid lines while those for modified measure are in dashed lines. 
The SR plots are shown in red. The figures in the left column are for Case I while those in the right column are for 
Case III. The DSR plots for Case I are plotted in green. For Case III, the leading order behaviours are plotted in blue. 
The numerical plots for Case III are also shown in green for comparison. Different values for the parameters in natural 
units are as follows: V = 10 35 , N = 10 25 , k = 1, m = 1 (Case I), 0.01 (Case III), a = 1, a% = 0.2, a 2 = a 3 = ... = 0. 
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In Fig 2, we have plotted the various thermodynamic quantities for Cases I and III as a function of T and compared 
them with the SR case. The P vs p plots have been obtained by varying T keeping all other parameters fixed. The 
qualitative natures of the plots for different cases are same. In case of S there are two competing effects: while the 
cutoff tries to reduce S by limiting the number of accessible states, the modified dispersion tries to increase S by 
enhancing the Boltzmann weight exp(— e/T) (note that £dsr(p) < £sr(p) for a given momentum state p, the change 
being more for greater value of the parameter m). At low temperatures, the latter is dominant and Sdsr > Ssr- For 
our choice of parameters this is clearly visible for the plot of S in Case I. In the high T regime, the cutoff effect comes 
into play and Sdsr < Ssr- The cutoff also saturates U as T increases, and C'v — > 0, resulting in a steeper equation of 
state. Here we make an interesting observation. There have been attempts to define velocity in DSR [34] . If we adopt 

the usual definition for the speed of sound c s = \J~^;, then we observe that c s grows without any bound. Possibility 

of such scenarios has been discussed in [35] . 

For given choice of parameters in case of quantities like S and U, the leading order behaviours for Case III match 
with the numerical plots up to T ~ 0.15, while in case of Cy which contains second order derivatives of the partition 
function with respect to T, the leading order behaviours match with the numerical plots up to T ~ 0.09. Note that the 
leading order behaviours have been obtained assuming -4— to be finite and ^ — >• which in turn implies 4- = — — > 0. 
Hence as T increases, the leading order plots depart from their numerical counterparts. 

We conclude with a summary of our investigation of the ideal gas thermodynamics in DSR framework. In this paper 
we have used the dispersion relation (1.1) and have considered the modified phase space measure (the modification 
being isotropic and expandable in Taylor series). We consider three cases separately (m = K,m > K,m < k). The 
single particle partition function has been shown to be smooth in m n € (0, n) (see Appendix B). For the case m = n, 
a simple analytical form for the partition function is obtained (see (2.17)) while a series solution for the partition 
function has been obtained for m < n (see (2.27) and (2.36)). In doing so, new type of special functions (Incomplete 
Modified Bessel functions) emerge. We observe that DSR thermodynamics is non-perturbative in the SR and massless 
limits. Using the leading order solutions, we derive thermodynamic quantities like the free energy, pressure, entropy, 
internal energy and heat capacity (see Fig. 2). A stiffer equation of state is found. 
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Appendices 

A Convergence of K* (x, y) 

K* (x, y) given by (2.36,2.38,2.39,2.41) is convergent if the following two series are convergent: 
and 

oo 2r-3 (-x) k ^ 1 fl\ 2r ~ 2 ~ k 

S ' = ^ 2 U g (2r-3)(2r-4)...(2r-2-fc) [yj " (A " 2) 

5 can be easily proved to be absolutely convergent using Ratio test. For S' first consider the following double series: 

oo oo / x ^k— 1 2r — 2 — fe oo oo 

S " = S (2r - 3)(2r - 4)...(2r - 2 - *) [y) ^S^' ^ 
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Let us first test the convergence of S" (See theorem (2.7) of [36]). The row series S r (for a fixed r) and the column 
series Sk (for a fixed k) are defined as 



oo 

S r ^ ^ &r,k? 
fc=l 

oo 

r=2 



The ratio tests for S r and Sk show that they are absolutely convergent (for y > 1). Also lim 

r,fc— >oo 

Hence, S" is absolutely convergent. Now 

oo oo oo 2r— 3 

|s"l = EEl fl ^l = EE IM + £ ;i>o. 



(A.4) 
(A.5) 
= < 1. 

(A.6) 



r=2 fc=l 



r=2 fc=l 



oo 2r — 3 

As 1 5" | is convergent (or in other words S" is absolutely convergent) we must have \ S'\ — ^ ^ \a r ,k\ to be convergent 

r=2 fe=l 

(or in other words S" to be absolutely convergent). Thus the series expansion of K* (x 7 y) is absolutely convergent. 

B Continuity and differentiability of the partition function in tuq 

We shall show that Z\ is continuous in too for m <G [0, n]. After integrating over the angular coordinates (2.13) gives 

Z " ( m o) - --^2 / d P p 2 e-^(p^)-m a] = v / dp /(pjmo) . (b.i) 

(2-7T) Jo (2-7T) Jo 



The integrand f(p, mo) = p 2 e -^[ £ (p. m o)- m o] is a continuous bounded function of p and mo in the range mo € [0, k],p£ 
[0, k]. Thus (mo) is a continuous function of mo as the function g(p) = C, where C is the upper bound of \ f{p, mo)\ 
in the range mo G [0,«],p e [0,/t], satisfies g(p) > \f(p,mo)\ for all mo G [0,/t],p G [0,/t] and is integrable as 

/ dp g(p) — Ck < oo (see Lemma 1 in § 5.12 of [32]). 

JO 

The derivative of the integrand with respect to m is given by 



df(p,m ) 
9m 



Pf(p, m ) 



e 2 (p,m ) -p 2 



1/2 



(B.2) 



It has 2 poles (and also branch points) in the complex p-plane at p = ± - 



, ... We note that the poles and the 

(^-i)' l 

branch points remain to be at the same positions for all higher order derivatives of f(p, m ) with respect to m . For 
m = both the poles are at p — and as mo increases the poles separate towards the imaginary axis. They keep 
on moving on the imaginary axis till they reach ±ioo at mo = After that they start to come closer to each other 
on the real line and finally at mo = k they stop at p = ±n. Note that for all mo G (0, k) the poles are never on the 

contour of integration (the real line from p = to p = k) and the functions ^ 8 ^q^"'°^ , 9 , etc.') remain to be 



bounded. This (by the same argument as given in the case of the continuity of Z\ (too)) ensures the infinite-order 
differentiability of Z\ (mo) in mo G (0, k) and the derivatives can be found by using the Leibniz rule (see Lemma 2 in 
§ 5.12 of [32]). Note that the fact that we are not being able to say about the differentiability of Z\ at to = could 
be a relic of the non-analytic part in (2.45). 
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